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Abstract 


This talk will first make an expository survey on the recent develop- 
ment of the low-Reynolds-number hydrodynamics pertaining to micro- 
organism locomotion. The scope of interest will be focused on two primary 
modes of swimming in this flow regime, namely the flagellar and ciliary 
propulsion. In order to overcome the shortcoming of classical resistive- 
force theory, the improved slender-body theory requires an accurate | 
account of (1) the long-range force field generated by the flagellar motion 
(which is invariably of finite amplitude), (2) the interaction between the cell 
body and flagella, and (3) the state of self-propulsion. These basic features 
interrelate to make the problem interesting, the development of the refined 
theory challenging, and the improved accuracy being achieved for theoretical 
prediction doubly rewarding to both the fluiddynamicist and microbiologist 
for their future joint studies. 


INTRODUCTION 


In nature the phenomena of locomotion of aquatic animals in aqueous 
media embraces a great variety of species with different sizes, different 
modes of swimming motion, and with varying biophysical basis of physiology 
and metabolism. Of these families of swimming organisms, the size spans 
over an impressive nineteen orders in magnitude; their basic mechanism of 
swimming ranges from the viscosity-dominant, inertialess fluid reactions 
for micro-organisms in one extreme to the inertia-based reactive hydro- 
mechanics for larger animals inthe other. In physiology, the metabolic 
processes vary from the better known conversion of chemical energy to 
mechanical work by means of muscular contractile effort in higher animals 
to the not-so-fully understood processes in prokaryotic flagella of bacteria 
in which no enzymatic activities are known to exist. Evidently this subject 
is interdisciplinary in nature, and its development would be most rewarding 
if cultivated by a joint team of bioiogists and fluid-dynamicists. 


The world cf micro-organisms that can perform swimming motion may 
be classified, according to their modes of movement, into three major 


categories. 


1. Flagellates. Flagellates, such as the various kinds of spermatozoa and 
many flagellated bacteria, locomote themselves by means of undulatory 
movement of a single flagellum or several bundles of flagella, in the forma- 
tion of planar waves (e.g. Ceratium, Ochromonas, Prymnesium), helical 
waves ‘e.g. Trichomonas, Trachelomonas), or of three-dimensional waves 
of a more general type (e.g. in Euglena and Peranema). Further variations 


of these basic forms include the spirilla, spirochaetes, and the flagellar- 
mastigoneme propulsion system. 


2. Ciliates. The ciliates are a class of protozoa which have a large num- 
ber of hair-like organelles, called cilia, attached to the cell surface ina 
fairly regular row-and-column distribution, performing a beat pattern in- 
herent to each species to form a metachronal wave in propulsion. Some of 
the common ciliates, whose motion has been more extensively studied, in- 
clude Paramecium, Tetrahymena, Spirostomum, and Opalina. The ciliary 
mode of prooulsion is found not only in ciliates, but also in some organs of 
higher animals such as the muco-ciliary systems in respiratory trachea and 
in reproduciive ducts. It is in the latter problem areas that further develop- 
ment of ciliary dynamics and physiology will have a most active future. 


3. Pseudopods. The temporary protrusion and retractile process of the 
protoplasm of a cell for moving about, such as in amoebic movement involves 
certain mechanism of intracellular flow and cytoplasmic streaming. This 
mode of movement appears to be the least understood of the three, and, for 
this reason, the present discussion will be devoted to the first two types of 
propulsive movements. 


From the hydromechanical point of view, these movements have three 
important features which are essential to our basic understanding. (I) First, 
the motion is generally characterized by two small Reynolds numbers, one 
based on the mean propulsion velocity and the other on the oscillatory compo- 
nent of the motion; they may be expressed as 


Rey = Ut/yv, Re. = ub? /y (1) 


where U denotes the mean propulsion velocity of a micro-organism of 
length £, uw the radian frequency of an undulating flagellum or cilium of 
cross-sectional radius b, and v is the kinematic viscosity coefficient of an 
aqueous medium (v = 0.01 cm2/sec.). Typically, Rey ranges from 10-6 
for small bacteria to 10-2 for most spermatozoa, and Re_ is approximately 
an order of magnitude smaller. At such low Reynolds numbers, the fluid 
motions so produced are dominated by the viscous force while the inertial 
effects, of voth the fluid and the microscopic body, are entirely negligible. 
This means that the flow would altogether stop as soon as the propulsive 
motion is halted, (II) Further, the velocity field generated by the motion of 
a body element has the longest range known in fluid mechanics and the motion 
is invariably of a finite amplitude. These features combine to make the 
hydromechanical analysis both interesting and challenging. The ratio of the 
amplitude, h, of a flagellar or ciliary wave to the wavelength, A, is typically 
about 1/7. The recent efforts have been directed primarily to improve the 
theory so that the effect of flageilar curvature can be more accurately deter- 
mined for motions of such large amplitude. (Il) In a broad sense, move- 
ment of prekaryotic flagellates and of eukaryotic ciliates and flagellates can 
be regarded as self-propelled, i.e. without aid of any extraneous force and 
moment of force. The mathematical expression of the state of self-propulsion 
is, symbolically, 

= forces = 0, and x moments = 0, (2) 
where the s‘gn £ includes the summation by the limiting process of integra- 
tion. This means that while we need to use forcelets and their moments to 
represent tle propulsive action of flagella and cilia, the resultant force and 
moment must always vanish, even when the micro-organism accelerates 
along a rectilinear or curvilinear path, because the inertia of both the fluid 
and the body is negligible. This implies that the fluid at some distance away 
from a self-propelling body will already feel ne trace of a net force and 
moment; it may sense at most certain higher moments of a force distribu- 
tion. Asa cesult, this feature brings forth some unique flow behaviors that 
are by their very nature absent in low-Reynolds-number hydromechanics of 
inanimated small objects. In combination, the above basic features of Now 
pertaining to micro-organism locomotion will form an overall basis for the 
yresent usscussion. 


ton micro-organism locomo- 
tion can be attributed timulation provided by the pioneering work of 
Sir James Gray (1928, 8) and Sir Geoffrey Taylor (1951, 19§2a,b). 

Recent advances in this active field have been reported extensively in the 
literature, e3>e y in the volumes by Sleigh (1973), Lighthill (1975a), Wu, 
Brokaw & Brennen (1975), and in the review articles by Jahn & Votta (1972), 


The rapid develo; 


and Brennen & Winet (197¢ 
FUNDAMENTAL FLOW SINGULARITIES 


In order to help us appreciate the recent development, in both its bio- 
physical and mathematical contents, of this general subject, it is useful to 
have an introductory disclosure of the basic tools which have been employed 
by the applied mathematician and hydrodynamicist for this class of problems. 


Insofar as the inertial effects are negligible, the hydrodynamics of low- 
Reynolds-number flows can be prescribed by the Stokes equations 


WVsw re, (3) 
Vp = py tiixt), (4) 


where u denotes the flow velocity, p the pressure, p the dynamic vis- 
cosit’ coefficient (1. = pv, p being the fluid density), and f represents a 
generalized distribution of extraneous force which may depend on the radius 
vector x andtime t. The primary fundamental solution of the linear 
system (3) and (4) corresponding to f being a singular point force in an un- 
bounded fluid was discussed by Oseen (1927), Burgers (1938), and called a 
Stokeslet by Hancock (1953). The higher-order poles associated with a 
Stokeslet have been explored by Batchelor (1970) Blake (1971-74), Chwang 
and Wu (1971-76) and others. A brief description of the properties of these 


flow singularities can be given as follows. 
(i) Stokeslet. — A Stokeslet of strength a(t) corresponds to having 


{ = 8mpa(t) d(x) , (5) 


—_ 


where 6(x) denotes the Dirac delta function. Its induced velocity and pres~- 
sure fields are, respectively, 


Uc(xia) = g/r+ila+x)x/r , (6) 

Polxia) = -24 0+ (a/r) , (r= {x{)- (7) 
It exerts on the fluid a net force 

Fs = Srp a(t) , (8) 


which requires an equal and opposite force acting on the fluid at infinity in 
order to maintain dynamic equilibrium. The velocity is seen to fall off like 
r-", a rate which is not surpassed by any other known flow singularities. 
The evidence of this long-range effect can be detected from the resulting 
streamlines (see Fig. 1), showing that fluid is drawn in from behind to 
approach the singular force and is then thrust forward along divergent paths 
with a fore-aft symmetry. 


(ii) Rotlet. — A rotlet of strength y(t) is associated with the singular 
moment of force os 


f= 4m VX yit) O(x) , (9) 


the resulting flow velocity and pressure being 


Fig. 1. The streamlines due toa 

Stoke slet of strength 2, shown ina oie agi es abe gs 
meridian plane about the axis of x 
symmetry. 


hee te to 2 rotlet 

4 visualized as the one d 
to a rotating sphere of small radius . 
with angular velocity L- 


Ug(x y) = OK (y/r) , Pp. = 0 (10) 


It exerts on the fluid a net moment 


Mp = 8mylt) , (11) 


which is balanced by an equal and o 

pposite moment on the fluid at 
Physically, the flow field of a rotlet may be visualized as the nak Fr 
small sphere of radius a rotating about the y-axis with angular velocity 


S 3 ; 
a= y/a (see Fig. 2). The streamlines are-axiall ymm 
the y-axis, with flow speed falling off like r-2. jig ae 


iii) St let. — 
ol “a et A stresslet is a component of the Stokes force-dipole de- 


Usp(e- B) =~ B+ WUe(eal = Un tUcs , (12) 


where Up represents a rotlet of strength yY = BXa, and 
‘ 3 
Usg(xia, B) = -(a- Bx/r> + 3a. xB- x) x/e?, (13) 


is the stresslet velocity field. This veloct 
. ty is symmetric with respect 
interchange of a and §, and represents a pure straining or ~buacienh tne 


(see Fig. 3). i 
= e ~ : oo the rotlet, the stresslet has its velocity decreasing like 


(iv) Mass 4oublet. — A mass doublet, of strength 6, has the same velocity 


field th i j 
porches at of a dipole known in potential Now; it is ‘also the Laplacian of a 


Up!iB) = V+ (B/r) 


Us(xip) . (14) 


By 
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stokes dipole rotlet stresslet 


Fig. 3. Decomposition of a Stokes dipole into a rotlet and a stresslet. 
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It has zero vorticity since Wp possesses a scalar potential and it carries 
no pressure on account of the negligible inertia effects. 


Higher-order flow singularities can be readily derived from the above 
fundamental cases; they are useful in construction of exact solutions (Chwang 
& Wu 1975) and especially when a high degree of accuracy is required 
(Chwang 1975). 


RESISTIVE-FORCE THEORY 


In the general case of flagellar movement, approximate calculation of 
the highly viscous flow past a slender body, such as a flagellum, can be 
carried out by applying resistive-force theory which seeks a definite rela- 
tionship between the force exerted by the body per unit length on the fluid 
and the local velocity of the flagella. The central concept of the theory is 
quite simple. Suppose a segment of a slender flagellum of arc length ds 
moves with velocity Vi{s,t) (relative to the unperturbed fluid far away) so 
that its component tangential to the body centerline is v and the comple- 
mentary normal component is V.. The resistive force éxerted by the 
segment on fluid has two components dF, and dF_ inthe same two direc- 
tions which can be expressed as a 


dF, = CypV,ds, aF = CinVids , 


where C_ and C_ are two dimensionless force coefficients assumed to be 
dependen? only on the instantaneous body geometry and independent of the 
segment position s orits velocity. For an elongated body of nearly circu- 
lar cross section, the force coefficients assume for a quite general category 
of motion the expression 


27 = 4n 
fe Tog(2q/b)-C, Cn* Toglaq/b) +c, ' (16) 


where b is the local or a characteristic cross-sectional radius of the 
flagellum end q denotes a length whose value depends, together with the 
two numerical constants Cc) and Cz, on body geometry as well as the type 
of flagellar movernent. For instance, for an elongated prolate spheroid with 
sermi-axes a ind b (b « a) (16) is accurate with a relative error of O(b/a)? 
(Tillett 1976) if 


q = a; c,=c = (17) 


The local force coefficients of a flagellum in undulatory movement are 
much more difficult to determine accurately. It was Hancock (1953) who 
first applied a line distribution of Stokeslet and mass doublet to represent a 
slender boay movement, and, based on a limiting result of his analysis, 
Hancock proposed the form (16) with 


q = 4/2, C,=C,= 5 (Hancock 1953) (18) 


1? 
where £ is the length of the filament. Subsequently, Gray & Hancock (1955) 
introduced for computing flagellar motions the formula (16) with 


(Gray-Hancock) (19) 


where \ is the wavelength of a planar or helical flagellar wave. The popu- 
lar appeal! of this formula is evidently related to the simplicity of the resulting 
mathematical analysis. The theory, however, is crude and its accuracy is 
insufficient for making more refined studies by our modern standard. Along 
this line of development Lighthill (1975b) has recently re-examined the hydro- 
dynamics of flagellar locomotion and suggested that for flagellar waves it is 
accurate to take 


26091, 2,8, &- 


: 2 => (Lighthill 1975p). (20) 

To provide an order estimate of the accuracy of these formulas, we re- 
call the theory for straight slender bodies (of length £4) developed by Tuck 
(1964), Tillett (1970), Batchelor (1970), and Cox (1970). Based on expansion 
of the line-force distribution in terms of a series in inverse powers of 
log (4/b), the theory shows that q actually depends on s for arbitrary body 
shape except for slender spheroids. But as an approximation, formula (16) 
with constant ae of Pages Cc, can still be used (e.g. q= 4/2, 
cy = +0. 807, Cc, Cds for a circular cylinder) if we accept that the error 
is of order [log” Us 709} With such errors, the accuracy of the crude 
resistance coefficients is generally regarded as insufficient for the scientific 
requirements and becomes even more questionable when the body centerline 
bends and waves as in flagellar motions. 


The need of having an improved slender-body theory was stressed and 
extensively discussed, among other basic problems, during the Pasadena 
Symposium on Swimming and Flying in Nature (Proceedings edited by Wu, 
Brokaw & Brennen 1975). This problem was further pursued at the 
Cambridge Symposium on Biodynamics of Animal Locomotion which was 
organized by Professor Sir James Lighthill and Professor Torkel Weis-Fogh 
(Proceedings edited by Lighthill & Pedley, to be published by Academic 
Press). From these rounds of discussion we may single out the following 
observations. 


Based on a series of experiments using mechanical models, Chwang and 
Wu (1975b) found considerable discrepancies between theory and experiment. 
Of the significant factors that may contribute to these discrepancies, we have 
subsequently identified: (i) the effect of finite curvature of body centerline, 
{ii) the wall effects due to the proximity of the container boundary (both in 
model and prototype tests), and (iii) the inertia effect that may arise when the 
Reynolds number based on the mechanical model length is no longer small 
(even when the body centerline is moving perpendicular to itself). 


New research has been pursued along these identified directions, The 
first test case is the slender circular ring (the only shape with a uniformly 
curved centerline yet without body ends). For the Stokes flow involving a 
ring in arbitrary motion we have found that the solution can be determined 
asymptotically to any specified order of accuracy, at least in principle, by 

ling higher-order flow singularities. This ring solution has been 
ructive in guiding us to establish a refined slender-body theory for the 


=> .-*® 
— _ ~~» 


\ \ f 


general case of arbitrary body motion with finite amplitude. Development of 
flagella and cilia hydromechanics, again no exception like many major fields 
of fluid mechanics, has been full of challenge, perplexity and excitement, 


Before I proceed with my report on this task, we conclude this section 
by noting that the force calculation must be cornplemented by the considera - 
tion of a pure couple (See Chwang & Wv 1971) acting on the fluid by a segment 


of flagellum whose surface undergoes a spin with angular velocity 2. paral- 
lel to the unit vector ge, directec along the body centerline. The corres- 
ponding moment is i 
2 
eM, = Cyn d Q,ds ; Cus = #2. (21) 


Unlike the force coefficients, this formule remains accurate as the body 
length is increased, whatever the value of the Reynolds number based on 
Q, and body length £ may be (Chwang & Wu 1974). 


THE STOKES FLOW PAST A SLENDER RING 
In order to investigate the body-centerline curvature effect in Stokes 


Now, Wu & Johnson (1976) have selected a slender ring as the primary test 
case. Itis of interest to summarize here the new result and to ascertain, 


by comparison, the deficiencies of the pr vious theories in predicting 
sectional forces and moments. The math«matica! details, however, will be 
curtailed since they are Leyord the pressn* scope 


The result is presented below separately for four basic modes of motion, 
namely (i) broadwise translation, (ii) edge-on translation ( in centerline 
plane), (iii) rotation about the central axis, and (iv) on-edge rotation. The 
general case of arbitrary motion can be resolved by superposition of these 
solutions. 


(i) Broadwise rotation — Consider a slender ring characterized by a small 
parameter « = b/a «1, where a,b are its major and minor axes, respec- 
tively (a being the radius of the centerline circle and b the cross-sectional 
radius). When the ring moves with velocity V parallel to its central axis 
(z-axis say), the following flow singularities, distributed uniformly along the 
centerline, are required to represent the flow, with their line densities: 


1 1 
Stokeslet (along He a = x VAL +>) = 
doublet (-¢,) B = i ba 
1 1 
rotlet (-¢9) y = a Vba(L-1)/(L + 5) (22) 
1 
stresslet (@,. e,) A= ba , 
1,2 
quadrupole {e.. ¢,) B= - q? A. 
Here, (e, fg) e,) denote the base vectors in the cylindrical polar coordi- 
nates (r, 6,7), and L stands for 


L* log(8/e) , « = bla«). (23) 
We note that the Stokeslet anddoutlet are required for the solution of flow 
velocity y upto O(L-") for any positive integer n, and inclusion of the 
next three singularities reduces the error of solution for yu to O(«Zlog ), 
which i> a considerable refinement. 


{ii) Edge on translation — 


Wh , , P 
Se bet spe en the ring moves in the x-direction (taken to 


¥ = &, * 5 * te" ¥. = Vcosé@, 4 = Vsiné, 
- (24) 
e required centerline distributions of leading singularities are found to be 
1 
Stokeslet (e.) o-F V(L-3)/(L2 - 3 L-l) , 
ms 1 5 
Stokeslet (-¢ ) a, = ZV,(L- sl? - 3-1), 
doublet B=- 7 a, (25) 
3 
rotlet (e.) y = gVgbee - Prey? Sn -1) . 
3 
stresslet (ee) A = 7 V,ba(L - By? - FL-1) , 
ad a 1,2 
guadrupole fee.) B= =3b A. 


The velocity field uy represented by th i i 
y these flow singularities is accurate 
-: @ small €: ror term of O (42 log «). However, it is of significance * note 
at in this case, the resulting force coefficients, C, = 8na,/V, 
a 


C_ = 8ra_/V no longer ass i 
= » ume the classical form (16); 
evidently due to the effect of the centerline ee es eee 


(iii) Rotation about the central axis — For the case when the ring rotates 


about the z-axis with angular velocity 9, we ha i 
distribution of Stokeslet and rotlet, with ‘ennai ihe ha cs 


Stokeslet (e @) a 


FMa/(L- 2), 


(26) 
rotlet (e_) Y 


3be aL - 4). 
In this case the error of y is also of Ole? log «). 


(iv) Edge-on rotation — In this mode the ring rotates about the x-axis with 


angular velocity Q so that at the instant when the centerl 
i 
the y= 0 plane the velocity of the centerline is ne 


Y¥*V se. Vv, = Qasine®. (27) 
The corresponding flow singularities have the following line densities: 
Stokeslet (e.) ae $V,ML - 3 ° 
doublet (-¢,) p = $b7e ; 
rotlet y = bea le, cot @ - HL - 2) £6) , (28) 
stresslet (e fy) A= 5 dea ‘ 
quadrupole (¢.,¢,) B= - 4? A. 


_ above results show that the sectional force coefficients have various 
orms cepending on the curvature and modes of motion. In applications to 


Magellar motion, log(f/b) usuall 
: sentiy rags: from 3 to 5 (corresponding to —_ 


i 


t/b = 20 to 150). For such moderate values of L, expansion of the 
rational! functions of L in the above expressions into series of L-! may 
induce unnecessary error. In this sense the present result can provide a 
useful guideline for developing the general slender-body theory to be dis- 
cussed below. Further, we note that the rotlet y will yield a pure couple; 
this contribution may be significant, especially when the net moment of 
Stokeslet force is small. The presence of rotlet may even cast new light on 
the mechanism of spirochete locomotion as discussed by Chwang, Winet & 


Wu (1974). 


To exhibit the difference between the new results and classical resistive- 
force theory, we present the numerical results in Figures 4-7 for the four 
basic modes of torus motion. The comparison shows that while the differ- 
ence is not large for the translational modes (less than about 15% for 
« < 0.1), it is quite significant in the two cases of rotation. This may be 
attributed to the feature that in rotation the body velocity has finite variations, 
both in magnitude and in direction, along the torus. On physical ground we 
may expect that changes in phase along a flagellum or cilium in undulatory 
motion would have a similar effect on the force coefficients. 
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piane. 


GENERAL SLENDER-BODY THEORY 


Johnson (1976) has considered the general case of a sl 
circular = i 2 = slender body with 
shaved id section b(s) of arbitrary distribution but with spheroidal- 


1/2 2 


o 2 4 2 
b= «(4° - 8°) for (1-s°/4£“) «1, e=bo/t«l, (29) 


where s is the arc length along the body centerline: 


x = x,(s, t) 


sa (-t < s< t), (30) 
the movement of x being assumed arbitra 
x ry. At each point 
centerline we attach a set of orthogonal base vectors (e nee oo ae 
re the tangential, normal, and binormal directions, respectively? (see Fig 
i.e. ‘ ; 


¢, ° dx,/as, &, ade ./ds » S&F e,Xe,: (31) 
where a(s,t) is the centerline radius of curvature, a = Jae /es|~? 
“es 


The velocity of material points on the body surfac 
7 e@isd - 
lation of centerline with velocity Ox,/at, a per spin preaepbesesti at 
with angular velocity Q(s,tle., Q. being arbitrary, and the rotation of 
centerline with angular Velocity ne = ex (8¢ ,/ a). Thus the no-slip bound- 


ary condition on the body surface is given by 


¥ 
a) - 


Kp + diel Oe, +8, X8,)Xs,_ (32) 


" 
< 
0 


ofe t Yntn * Vy Sy 


The geometry and coordinates adopted fora slender-body in 


Fig. 8. 
arbitrary movement. 


where the dot denotes differentiation with re spect to t, and 


cosy+e, siny - (33) 


e = Sp 


£n 
resented by the 


The solution of the resulting Stokes flow can agin e body centerline 


distribution of necessary fundamental singularities on th 
as follows: 


° . 
ulx.t) = { (Ug (Rig) + Up( BiB) + UplBiy) + Usq(Rim) + Ay Uss (Rie, &p) 


-s 
° 


J > = sels ds, 
~ BU (Riegs Sn) + AgUsg (Riggs Sp) - Baa (Bien e,)} 


(34) 


h 
sie R=x-x(s) . R= |R| , (35) 


U_(R;m) = mR/R° is the contri- 


»U,U,U have been given before, U_(R : ie : 
EMisPor BinSSs source of strength m, and U (R69) = (6 vie he 
denotes a quadrupole. The limits of integration, + S,, are g 
generalized foci, 

8 ey ~~ 508 - (36) 


°o 


i i hood of the body (i.e. for 
in (34) is expanded in the neighbor e 
Bais ck nach etation Pigg integration and then applying om grad 
Saadition (32), neglecting the terrns of Olae log «), we obtain the following 


integral equation for g, 


56 


° 

Vi (s, t) = a(s,t)L +f K (Roads! vs a; (37) 
-3 
°o 
where 
L, = 2[2 log(2/¢,)- 1) Li = Ly =2 log(2/«,) aie: 
a,(s',t}  (a(s', th RL)R Da (s,t) 
K (Bia) = == i - ; (38) 
sera R R |s-s'| 
° ° 
D, = 2, dD, = dD, a er &,= b(s)/t(1-s2)!/2 : 

R, = x, (8, t) — 2 x,(8', t) = Rosfs + Rontn + Ro» Sb x 


(39) 


Here, the base vectors e.,e., fy in the expression (37) assume their 
value: at s where V,(s,t) is specifier. The integral eq ution (37) can be 
solved by applying certain appropriate iteration schemes or the method of 
collocation, These numerical methods have been applied to evaluate sec- 
tional force coefficients of a flagellum performing a planar or a helical wave 
motion. A similar theory has been introduced by Keller & Rubinow (1976) 
for the leading singularities, namely the Stokeslet and doublet, but with their 


line densities only approximately determined near the body ends. 


WALL EFFECTS 


Low-Reynolds-number flow theory has applications in several fields of 
growing importance, for instance in biophysical and biomedical applications, 
rheology, suspension mechanics, aerosol physics, and micro-organism 
motions. In these applications it is often necessary to analyze the motion of 
slender bodies such as fibers or organisms in the vicinity of solid boundaries 
since their motion is considerably modified due to the wall effects. Of 
these phenomena some appear in natural scheme of life (e.g. the muco- 
ciliary propulsion system in trachea and spermatozoa transport in ciliated 
duct), some are results of artificial design (e. g. the flow of suspensions 
during engineering processing), while many others are caused by necessary 
preparation for microscopic observation. We thus see the importance of 
wall effects in low-Reynolds-number flow, both from the standpoint of basic 
understanding and data interpretation. The general problem, however, is 
difficult, and as a result, only a few relatively simple cases have been 
considered. 


The image system for a Stokeslet and higher-order poles into a plane 
wall can te derived from the classical work of Lorentz (1907), as illustra- 
ted by Burgers (1938), and more systematically identified by Blake (1971a, b, 
1974), Blake & Chwang (1974). This method has been applied by Brenner 
1962) to discuss arbitrary bodies, and by DeMestre (1973) for the case of 
circular cylinders near a plane wall, by DeMestre & Russe} (1975) to study 
rotating rod motions. Further applications have been made by Winet (1973), 
Blake & Sieigh (1974), Katz (1974), Katz & Blake (1975), Katz, Blake & 
Paveri-Fontana (1975), and by Aderogba (1976) for various studies of ciliary 
and micro-organism propulsion, A review of these contributions can be 
found in Brennen & Winet (1976), As an overall summary of these results, 


a - ~ eee Cece. ae a 


the following feature is noteworthy from both the hydromechanical and the 
energetics points of view. 


Generally speaking, the drag increases as the wall is approached. But 
the relative increase in drag of a slender body is always smaller in axial 
motion than that when the motion is transverse to the longitudinal body axis, 
whatever its inclination with respect to the nearby wall, ie. 


©, rn +6), oS C0! +6, 


o< é.< bo 

(40) 
where C and Cyo are the axial and normal force coefficients of a given 
body in ufBounded flow, (i.e. for vanishing ratio of the body length to the 
distance from a wall, £/h = 0), and the relative increments 6. and 8. 
depend on X = f/h and on the body orientation with respect to the wall.” The 
ratio of om to ome is, by (40), 


y= O/C, 2(C JC, IFS WIL +5) < O/C, = y,- 


(41) 


rally reduced by the wall effect, 
43 pointed out by Lighthill 

n important role in swimming 

ird velocity, V. and to the 


max, 
Triax’ * 


Thus, the ferce-coefficient rati 

and this is a result of especial sig 
(1969), the force-coefficient ratio y 
propulsion. It relates to the maxim 


flagellum according to 


o 


maximum efficiency of locomotion 


v / =(1- (42) 


c 
max b 


n 
‘max : 


where c, is the phase velocity of flagellar wave with respect to body, and 
is? dat V/c, = 1-/fy. Inunbounded flow, y = 1/2 is the 
Yinan °* reache at Vic ‘Y- ‘ ay: 
lower limit reached only by extremely slender spheroids, and, as pointed 
out by Wu (1975), no swimming in the microscopic world would be feasible 
if y could attain the value y= 1. in most cases of micro-organism loco- 
motion, y ranges between 0,6 and 0.7. In this rather narrow range of y 
in practical operation, we thus see (from(41) and (42)) the possible advantage 
for a slender body to swim near a wall since the value of y can thereby be 
somewhat decreased (by about 10 - 15%), giving both higher swimming veloc- 
ity and propulsive efficiency. 


The general problem of wall effect can be formulated for the case ofa 
plane wall by incorporating the known image systems. Suppose that the 
velocity of an unbounded Stokes flow due to a slender body with given move - 
ment is known, which can always be expressed as 

8 
° 
. iy {2g(xis,t,)- V[(x-¢)- g)+Ve pds, (43a) 


-s 
° 


where g = (91+ %2 93) and %, are harmonic functions and ¢ a constant 


= 


vector, 
v¢ = 0, v6, = 0, & = (C),Co, (43b) 


Then the flow velocity generated by the same body sha; ent in 
the presence of a wall coinciding with the x, = 0 plane : 


o 
B -f {2p (xis,t)- V[ie-o)- gy) + Wjds (x, 20), (44) 


$ 9 = P(x) - Ax) » G(X] + XQ %3) = wo ley, X2, - Ky), 


e > oth - Sia + 28 sla (- (x- 9) g)+ + (x,9)) . 


P(x) + X>, %3) = 9(X),%>, - =) ; 


In the integral representation (43) of the unbounded flow, it is easy to identify 
the different terms in the integrand of (43) with various flow singularities. 
The representation (44) is obtained by appropriate imaging so that y = 0 at 
the wall x,=0, and y behaves like Uy near the body. By expanding the 
integrals in (44) for the near field and applying the boundary condition (32), 
we can convert the resulting integral equation into a form similar to (38), 

and it can therefore be solved by the same iteration scheme. Numerical re- 
sults are being analyzed for a few typical cases. 


INERTIAL CFFECTS 


Inertial effects will arise to effect changes in sectional force coefficients 
when either of the Reynolds numbers defined in (1) becomes no longer small 
compared to unity. Such situations may occur in motions of large micro- 
organisms or in experiments with mechanical models whose corresponding 
Reynolds numbers cannot all be kept small in practice. 


A clear explanation of how inertial effects can arise has been provided 
by Chwang & Wu (1976) for a prolate spheroid, Even when the translational 
motion is directed only along the minor axis, with velocity U say, the 
solution is found to depend nevertheless on two Reynolds numbers, F_=Ua/v 
and R, = Ub/v, based on the major and minor semi-axes a and b, Fespec- 
tively. The solution based on an approximate singular-perturbation method 
is shown to agree with the formula of Oberbeck (1876) for the Stokes flow as 
long as both R, and are small compared to unity, but deviates from 
itas R, becomes of order unity, and finally approaches the Oseen's 
solution for a two-dimensional cylinder as R_ +o, being kept small 
throughout. This gradual transition thus provides a clear physical picture 
and explanation for the manifestation of the 'Stokes paradox' known in 
viscous flow theory. 


It is thus important to have such inertial effects accounted for in making 
theoretical predictions and in interpreting experimental results. Develop- 


ment of the general slender-body theory along this line is still in its early 
stage, but is noted to be underway. 


Recent development in cilia propulsion theory has been based on the 


‘envelope model’, the 'sublayer-mo ', and the 'traction-layer model'. The 
historical background and principal lts have been extensively dis- 
cussed and reviewed by Blake & Sleigh (1474), Keller, Wu & Brennen (1975), 
Wu (1976), Brennen & Winet (1977), so 1 shall not repeat here what is al- 
ready available in the literature. However, I would like to apply again the 


"traction-layer rnodel' to make performance predictions on the basis of 
optimization calculation. 


The traction-layer model is based on the concept that the discrete 
forces of a system of cilia can be represented by an ‘equivalent’ continuum 
distribution of an unsteady body force within the volume of the ciliary layer, 
provided the cilia distribution is sufficiently dense, Conversion from the 
discrete to the continuum force is based on resistive theory which relates 
the original sectional force of a cilium to its velocity relative to the flow 
produced in its neighborhood by the entire ciliary system. 


Using this procedure the body force within a plane ciliary layer can be 
expressed as 


in{kx-wt) 


f(x, y,t) = (O<y<t), (45) 


TI 


N 
= five 


° 


where N < wo, = 27/k is the metachronal wavelength, w is the radian 
frequency, and f vanishes at the cell surface (y = 0) and also outside the 
cilia layer (y >£), 4£ being the cilia length. The distribution functions 
f(y) can be determined by carrying out the conversion from the discrete to 


continuum force, Since f (y) are the Fourier coefficients of the continuum 


force at fixed y, the above expression for f holds valid for the general 


case of ciliary movements of finite amplitude. 


The leading term of the expansion, f (y), denotes the mean force distri- 
bution; it plays the essential role of providing the mean velocity field, 
uly) which satisfies the equation 


~ 
ay 
In fact, the mean propulsion speed of the cilia layer is given by the value of 
u_ for y>#. The two functions {| and u_ also determine the mean rate 


of working by fy on the fluid in mafntaining the mean velocity field, which is 


t 
E -f filyuyly) dy - (47) 


f(y) - (46) 


° 
° 
On the other hand, the rate of working by the cilia in producing fo is 
L 
en { f(yacty) dy» 

° 
where uly) is the x-component of the mean velocity of cilia. Based on 
resistive-force theory, the ciliary force is proportional to the flow velocity 
relative to the ciliary segment, i.e. f (y) = »C(u_-u_), C being a force 


coefficient. It therefore follows that the Gf extremizing Ey for 


— E, for fixed mean square of {fy or 


i 
8 
p, =f {o (y) dy . (48) 
° 


Now the first variation of E 


gives a= AD.» \ being an undetermined multiplier, 


1 
5(E, - XD.) = if {(u, - 2X.) bf (nm) + £4(n) bus(n)} dy =0 
° 


where m= y/£. Substituting (46), or 6f = -2,2 2 
: ‘ z apes d'(6u_)/d in the ab 
integral and integrating by parts, we obtSin an hate Nola EA — 
, f 
8(E, v AD.) = 0 the Euler equation e first variation 


y+ eee (Og nel), (49) 


Su, = 0 and (du, /dn) =0 (at mn = 0,1). (50) 


i 
ples ing £,(0) = £0) = 0 is the eigen- 


f(a) =pAsinksg , k= (tual! soe, (51) 


and th i 
n e relevant solution corresponds to n = 1 k, =. The resulting 


: 
sof ws yee 8 (46) jointly with (51) under condition u,(0) = 0 


uo(n) = AL Ia)? (™ + sinnn) (O<H<1l), (52) 


and (n) = ae D = . 
2 wu, yet? sant u() for n>. If Ds = ~ = 1 by normalization, 
en A =/2, and the corresponding E, = Ej/ut” = (Alm)? = 3/e* s 0. 304. 
The ratic E /B = pe /D #* = 0, 304 i 
; = = 0. gives the minimum energy impart 
to fluid for fixed? mean°squa re force; this value is smaller than ioe neanoineatl 
Ne eign: of the optimum solution obtained previously by Wu (1976) (in the 
atter case rele = 0. 364, corresponding to about 20% more of the mini- 
mum specific work done), The mean force and mean velocity profile, 
eer a * - 
{,/f, = 2* sinnn, uo/4, * 24, 289 +sintn), (53) 
are shown in Fig. 9. (Here, in phy i i a 2 
iz : . physical dimensions, f 
Uy = w/k is the metachronal wave velocity.) The slope°of er ne 


obtained by i : 
itaxtel by integrating (46) across the cilia layer, is related to the force 


dy z t {ly)dy , (54) 
° 


which asserts that the viscous skin friction at cell surface is counterbal- 
ig the fluid force acting on the cilia over the same area of cell surface. 
The higher-order forcing functions f{ (y) in the expansion (45) can also be 


optimized in a similar manner, but we shall not go into this here in 


preference 


of three-dim 


0.8 
y/, 
0.6 
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Fig. 9. The optimum distribution of the average ciliary-continuum-force, 


{_, giving minimum power required for fixed mean square force, and the 
corresponding velocity profile. 


STATE OF SELF-PROPULSION 


We have just seen how the ciliary force can be optimally distributed : 
when the curvature of the cell surface is zero. To extend this investigation 
to bodies of three-dimensional shapes, Keller & Wu (1977) have arg a 
simple hydrodynamic model to determine the self-propulsion of otiiaye 
micro-organisms and to examine the effect of body shape upon =e scary a 
of propulsion. is model is a further generalization of the envelope mode 
adopted by Blake (1971) and Brennen (1974) for considering the sah = 
of spherically shaped micro-organisms (originally based on the seh pee 
that the fluid adheres to a material envelope of the ata sbiomtosicciaer and is 
closely analogous to the control-surface model of Blake (197 3). ——- , 
present control-surface model we c ynsider the class of cilia ape ba 
prolate spheroidal shape, and assume that the mean flow velocity may have 


+ 


’ Ar 1 ee eee waRe, 
at the cilia envelope both a normal and a tangential compo 


uw. = V_ sid wu. ==V_cos6, (55) 
38 8 : n n 
i he di a 
where V are two constants and @ is the angle Ustween = direc 
S and the unit outward normal n to the spheroidal control 


tion of p ie 
, ; 10). The above simple relations represent the leading 

eneral variation in u_, and u long the cilia envelope; 

f } 1vior (as to have a normal flow velocity moving 

; > 


vation of Cheun 
have been supported by the observation of Cheung 


Fig. lo. The rate of enesgy cost for a prolate spheroid of constant volume 
and different eccentricity e. On the left the geometry of the control- 


surface model is shown with the cell body, the cilia layer and the control 
surface S. 


The Stokes flow satisfying boundary condition (55) can be constructed by 
an axial aistribution of Stokeslet and mass-doublet. The resulting energy 
expenditure, E, is shown in Fig. 10 for spheroids of varying eccentricity e, 
the results being given wth respect to the reference value of ér.U2a, 
where a and b=a(l-e yl/2 are the major and minor semi-axes, respec- 
tively. It is of significance to notice the following effect of body shape on 
energy expenditure. For the spherical shape (e = 0), the rate of working in 
moving a rigid sphere, of radius a, at velocity U is at the cost of an 
extraneows force F sothat E = FU = 67,.U2a. In comparison, a ciliated 
sphere with zero normal component of flow velocity (Vq = 0) and with the 
tangential component (V,) so adjusted as to have zero force (with the net 
Stokeslet strength a = 6) must spend energy E = 12p,U2a in maintaining 
the flow, which is then irrotational. Although the energy expenditure in the 
latter case is twice that for the rigid sphere, it does make a drastic distinc- 
tion that the sphere with the slip-flow is capable of propelling itself! When 
the spheroid becomes slender, the energy required for propelling a spheroid 
with a slip flow actually turns out to be less than that for translating a no- 
slip spheroid of the same eccentricity e for e > 0.82. With the additional 
freedom of suction and blowing (V_ # 0) the energy cost can be further re- 
duced (see, e.g. the curve with vA/U = 0.2 in Fig. 10). 


Streamlines of the resulting velocity field in the lab-frame are shown 
\) for the two basically different cases, one with slip-flow and the 
isfying the no-slip condition. In Fig. 12 two time exposures of two 
cia, one self-propelling whereas the other being immobilized and 


ie) 


The streamlines (in the 
y frame) of a self-propelling 
centricity 0.9 


Fig. 12. Streak photographs of a 


freely swimming specimen of 


(upper) and 


Paramecium caudatum 


of a dead of 


ing (towards 


right) under 


ress.on 


n the point I wish to make, namely; while there wil 
always be stimulating problems along the frontier where the scheme of life 
is carrying or ne must be careful in scientific experimentation so as not 
sment of life in its true manifestation. 


to lose the basic ¢ 
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